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Abstract 

Weak-coupling expansions (conserving approximations) are carried out for 

the attractive Holstein and Hubbard models (on an infinite-dimensional hy- 

percubic lattice) that include all bandstructure and vertex correction effects. 

Quantum fluctuations are found to renormalize transition temperatures by 

factors of order unity, but may be incorporated into the superconducting 

channel of Migdal-Eliashberg theory by renormalizing the phonon frequency 

and the interaction strength. 
Pacs:74.20.-z, 71.27.+a, and 71.38.+i 
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Interacting electron-phonon systems sustain both charge-density wave (CDW) and su- 
perconducting (SC) order. Migdal and Eliashberg || showed that a self-consistent theory 
for the electron-phonon problem can be constructed that neglects so-called vertex correc- 
tions (Migdal's theorem) in the limit were the phonon frequency is small compared to the 
fermi energy. Numerical calculations of SC transition temperatures based upon experi- 
mentally extracted (interacting) electron-phonon spectral densities [|J soon followed. The 
Migdal-Eliashberg (ME) theory is successful for predicting transition temperatures (and 
other materials properties) of virtually all low temperature superconductors @||. 

However, ME theory assumes that the electronic bandwidth is infinite, and that vertex 
corrections may be completely ignored. Although the former can be handled by simple 
modifications of the basic theory the latter require complicated extensions of the ME 
framework. Grabowski and Sham |7| carried out such an analysis for the interacting electron 
gas and found that vertex corrections strongly suppressed T c for plasma frequencies larger 
than ten percent of the fermi energy. Van Dongen || analyzed the Hubbard model (which 
can be viewed as the infinite-frequency limit of an electron-phonon model) and showed that 
quantum fluctuations (vertex corrections) reduced Hartree-Fock transition temperatures by 
more than a factor of three. 

In this contribution, ME theory is extended to include vertex corrections and produce 
weak-coupling expansions for T c that are accurate for the entire range of phonon frequencies 
(Coulomb repulsion effects are neglected here). Comparison of these weak-coupling expan- 
sions are made to quantum Monte Carlo simulations over a wide range of parameter space. 
Effects due to vertex corrections are strong when expressed in terms of the bare microscopic 
parameters, but can be incorporated into the ME theory (to lowest order) by renormalizing 
these parameters. 

The electron-phonon Hamiltonian chosen here is that of the Holstein model || in which 
conduction electrons interact with localized (Einstein) phonons: 

H = ~rn £ (SW + cl a c ja ) + </5>;(n iT + n n ) + \mU? + ^£p| (1) 
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where c] CT (cj a ) creates (destroys) an electron at site j with spin a, rij a = c^ a Cj a is the electron 
number operator, and Xj (pj) is the phonon coordinate (momentum) at site j. The hopping 
matrix elements connect the nearest neighbors of a hypercubic lattice in ^-dimensions and 



is 



for concreteness the infinite-dimensional limit (d — > oo) of Metzner and Vollhardt [10 
taken (1/d corrections should be less than 5% in three dimensions ||). The unit of energy 
is chosen to be the rescaled matrix element t*. The phonon has a mass M (chosen to be 
M — 1), a frequency Q, and a spring constant k = MQ 2 associated with it. The electron- 
phonon coupling constant is denoted by g; the effective electron-electron interaction strength 
is then the bipolaron binding energy 

U = --^-- 9 - (2) 

~ MQ 2 ~ k { ) 

In the limit where k remains finite and Q is large compared to the bandwidth (O — > oo, U = 



finite), the Holstein model maps onto the attractive Hubbard model [II 



t* 

H = -77^ ( c h Ck ° + c L c i-) + U H n 3\ n 3i ( 3 ) 

with U defined by Eq. (g). 

Our weak-coupling expansions are based upon the conserving approximations of Baym 
and Kadanoff []T^ ]: the free energy $ is expanded to a given order in the interaction strength; 
the self energy S(zu; n ) is determined by functional differentiation E(iu; n ) = 8§ / 5G{iuj n ) at 
each Matsubara frequency u n = (2n + l)nT] and the irreducible vertex functions T(iu m , iu n ) 
(in the relevant channels) are determined by a second functional differentiation. The free 
energy must be expanded to order U 2 in order to determine the correct transition temper- 
ature in the limit \U\ — > ||. The necessary diagrams for the free energy in the Holstein 
model are presented in Fig. |]. Note that conventional ME theory is identical to a first-order 
(Hartree-Fock) conserving approximation for the SC channel (with the exception that the 
phonon propagator is dressed), but is quite different in the CDW channel, because exchange 



diagrams are important even at the Hartree-Fock level [12 



The many-body problem in infinite dimensions maps onto a self-consistently embedded 



Anderson impurity problem because the self-energy has no momentum dependence 1 



The self-consistent embedding of the impurity problem is solved by iteration [[14]]: the self 



energy is calculated from the dressed (local) Green's functions by a conserving approximation 
for the impurity problem and the new (local) Green's function is determined by the 
embedding 

G(icu n ) = r dy , P ^ EE F^iUn + Z n [G}) (4) 

J-oo tUJ n + /! — 2j n [Lr \ — y 

with p(y) ee e-xp(—y 2 ) / tJii the bare density of states in infinite dimensions. 

To illustrate the connection between the microscopic parameters of the Holstein model 
and the renormalized parameters of the ME theory, the SC transition temperature is calcu- 
lated in a square-well approximation at half-filling. The transition temperature is assumed 
to satisfy T c « min(f2,t*). In this case, the electron self energy becomes (to lowest order 
in U) [0 

ee iu) n Z{iu n ) ; Z(iu n ) = l-U dyp(y) — - (5) 

Jo u* + (Sl + y) 2 

and the irreducible vertex function in the SC channel satisfies (to second-order in U) |12| 



"sc(iu m ,iu n ) = UT{9(u c - \u m \)9(u c - |^„|)[1 H / dyF^(iy) 

n Jo 



n 2 + y 2 



/ dyF^(iy) — } (6) 

7T jo (il z + y z ) z 

in the square-well approximation [|J with uj c the cutoff frequency. The transition temperature 
is determined in the standard fashion ||-||] 

~>2 / • m r o r j p(y) n i 



Tc/r = eM 'Wm ]xexp[ ~7^i ^^)] xex p[- 2 y % ( o)(fi + ,) 2 j 2 

x{0.85exp[-^ r^^M tan_1 -]> Xex P^ r dyF ^0&^^ 1+1 9 0&-^ W 
n Jo y p(0) uj c V7r jo \¥ + y z 2\V + y z 



including all nonvanishing factors in the limit \U\ — > 0. The first factor is the contribution 
to lowest order in \U\. The remaining terms are constant factors arising from the following 
four sources (in order of appearance): the phonon self-energy contributions; the electronic 
self-energy contributions; the finite bandwidth contributions; and the vertex corrections. 
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The phonon self-energy term (equal to exp [2.493] in infinite dimensions) is normally incor- 
porated into the dressed phonon propagator in ME theory A = |C/|p(0)/[l — 2.493|£/|p(0)]. 
The electron self-energy term approaches e _1 as Q — > and e° as Q — > oo. The finite- 
bandwidth term should be proportional to Q if the phonon frequency is much smaller than 
the electronic energy scale, but should be replaced by the electronic energy scale when the 
phonon frequency becomes the largest energy scale. For concreteness, the cutoff frequency 
is chosen to be uo c = 0.6Q, so that this factor approaches the correct limits of 0.69Q as 
Q — > DH and 0.85 as Q — > oo [|J. Finally, the fifth term contains the vertex corrections 
which vanish as O — ^ and approach exp [—3.739] as Q —>■ oo ||. The above equa- 
tion for T c has the correct limits as Q — > [T c = 0.69fiexp{ — (1 + A)/A}] and as Q — > oo 
[T c = 0.24exp{-l/|[/|p(0)}]. The coefficient of exp[-l/|l7|p(0)] in Eq. (g) is plotted in 
Fig. [| Note that the maximum occurs at Q ~ 0.8t* (approximately 20% of the effective 
bandwidth), but the transition temperature is weakly dependent upon Q for larger values. 
The ME form for the transition temperature, T c = 1.14u; c exp[— (1 + A)/A], can be fit to 
the above form in Eq. (0) by renormalizing the electron-phonon interaction A and the cut- 
off frequency uo c (since the contributions from second-order terms [phonon self energy and 
vertex corrections] simply modify the prefactor of the T c formula), explaining its robustness 
for real materials calculations. 

In order to make comparisons with the quantum Monte Carlo calculations [0, the 
self-consistent perturbation theory is solved numerically. An energy cutoff employing 256 
positive Matsubara frequencies is used, and the self-consistent equations are iterated until 
the largest deviation of each G(iu n ) is less than one part in 10 8 . At half-filling the instability 
always lies in the CDW channel. The effective electronic bandwidth in infinite dimensions 
is of order 4t* and the phonon frequency is chosen to be approximately one-eighth of this 
electronic bandwidth (Q = 0.5t*) as in the quantum Monte Carlo simulations. 

Both the first-order and second-order conserving approximations are compared to the 
quantum Monte Carlo data in Fig. ^. The transition temperature determined in the second- 
order approximation is found to be renormalized (by more than a factor of two at low 
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temperatures) from the transition temperature calculated in the first-order approximation. 



The second-order conserving approximation deviates [TIJ from the Monte Carlo results when 
g is greater than OAt*. Furthermore, the turnover of the transition temperature curve as 
g increases is not reproduced by these finite-order conserving approximations. Probably an 
infinite-order summation scheme (such as the fluctuation-exchange approximation JT2J) will 
be required to properly account for this turnover. 

In the high-frequency limit (Q — > oo), the Holstein model maps onto an attractive 
Hubbard model. In this limit many "direct" diagrams cancel against "exchange" diagrams 
(there is no interaction between electrons with the same spin) and the perturbation theory 
is easier to control. Quantum Monte Carlo simulations have already been performed ||14|| , 



and weak-coupling expansions through second order (for transition temperatures) have also 
been investigated with different techniques PJL7|. 



The relevant diagrams for the free energy of the Hubbard model through third order are 
depicted in Fig. The transition temperature in the CDW channel is plotted in Fig. |5| for 
the first through third-order approximations and for the quantum Monte Carlo simulations. 
Note that the agreement with the Monte Carlo data is excellent and that the weak-coupling 
expansions appear to have an alternating character for the Hubbard model with odd orders 
overestimating T c and even orders underestimating T c . 

Weak-coupling conserving approximations (that are valid for all values of the phonon 
frequency) have been carried out for the Holstein (and Hubbard) model at half-filling. Cal- 
culations must be performed to second order in the effective electron-electron interaction 
in order to produce the correct limiting behavior as U — > 0. Comparison was made with 
quantum Monte Carlo simulations and excellent agreement was found for weak coupling. 

Migdal's theorem states that, in the small-frequency limit, transition temperatures may 
be determined with a first-order approximation (using a dressed phonon propagator) because 
the vertex corrections vary as Q/t*. This result is easily seen in our framework, but the 
corrections to T c due to a nonvanishing phonon frequency can be large (when expressed in 
terms of the bare parameters). Migdal-Eliashberg theory is, however, quite robust, and the 



effects of a finite band-structure and vertex corrections (to lowest order) can be incorporated 
by renormalizing the phonon frequency and the interaction strength (in the superconducting 
channel) . 

Extensions of these results off of half-filling, comparison with other techniques such as 
the fluctuation-exchange approximation |T2|] or the iterated perturbation theory |TB[, and 
an investigation of Coulomb effects will be presented elsewhere. 

We would like to thank E. Nicol for many useful discussions. This research was supported 
in part by the NSF under Grants No. PHY89-04035 and DMR90-02492. 
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FIGURES 

FIG. 1. First and second order contributions to the free energy in a conserving approximation 
for the Holstein model. The solid lines are dressed electronic Green's functions and the wiggly lines 
are phonon propagators. The spin indices a and r are summed over. 

FIG. 2. Constant factor in the analytic formula for the superconducting transition temperature 
in Eq. (7) plotted as a function of phonon frequency. Note the constant term has a maximum at 
n « 0.8t*. 

FIG. 3. Charge-density-wave transition temperature for the Holstein model at half-filling with 
ft = 0.5t* plotted against the electron-phonon interaction g. The first order (solid line), second 
order (dashed line), and quantum Monte Carlo results [15] (solid dots) are included. 

FIG. 4. First through third order contributions to the free energy in a conserving approxi- 
mation for the Hubbard model (the dashed line is the electron-electron interaction U). The spin 
index a is summed over. 

FIG. 5. Charge-density-wave (and superconducting) transition temperature for the Hubbard 
model at half-filling plotted against the interaction strength \ U\. The first order (solid line), second 
order (dashed line), third order (dotted line), and quantum Monte Carlo results [14] (solid dots) 
are included. 
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